Spin-selective Kondo insulator: Cooperation of ferromagnetism and Kondo effect 
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We propose the notion of spin-selective Kondo insulator, which provides a fundamental mechanism 
to describe the ferromagnetic phase of the Kondo lattice model with antiferromagnetic coupling. 
This unveils a remarkable feature of the ferromagnetic metallic phase: the majority-spin conduction 
electrons show metallic- while the minority-spin electrons show insulating-behavior. The resulting 
Kondo gap in the minority spin sector, which is due to the cooperation of ferromagnetism and partial 
Kondo screening, evidences a dynamically-induced commensurability for a combination of minority- 
spin electrons and parts of localized spins. Furthermore, this mechanism predicts a nontrivial 
relation between the macroscopic quantities such as electron magnetization, spin polarization and 
electron filling. 

PACS numbers: 71.10.Fd 71.27.+a 71.30.+h 75.20.Hr 



Even 30 years after their discovery heavy-fermion sys- 
tems attract much attention due to their fascinating 
properties. Apart from being Fermi liquids with ef- 
fective mass thousand times as large as the free elec- 
tron one, they show all kinds of competing or coexisting 
phases, and at the boundaries between these phases one 
frequently observes quantum phase transitions, accom- 
panied by barely understood non-Fermi liquid behavior 
[l|-|3|. Heavy fermion compounds usually include lan- 
thanides or actinides with open 4/- or 5/-shells, which in 
the simplest theoretical modeling can be viewed as a reg- 
ular lattice of local moments coupled to the conduction 
electrons. This coupling typically leads to two competing 
mechanisms: the long-ranged RKKY interaction and the 
local Kondo screening. While the RKKY interaction fa- 
vors a magnetically ordered state, the Kondo screening is 
usually considered to form a paramagnetic heavy-fermion 
state. The competition of these two mechanisms can be 
easily understood in terms of the Doniach phase diagram 

i- 

While in most heavy-fermion compounds the magnetic 
order is antiferromagnetic, there are a certain class of 
compounds showing ferromagnetic order. For example, 
the recently discovered YbNi4P2 is a ferromagnetically 
ordered heavy-fermion compound which seems to be very 
close to a quantum critical point [5]. Taking such fer- 
romagnetic heavy fermion compounds as motivation we 
analyze in detail the mechanism stabilizing the ferromag- 
netic state. An interesting question in this context is, if 
and howjhe Kondo effect accounts for the ferromagnetic 
state 



In this letter, we propose a spin- selective Kondo in- 
sulator, where the Kondo screening plays an essential 
role in stabilizing the ferromagnetic metallic state at zero 
temperature, which elucidates a previously unrecognized 
feature of the ferromagnetic phase: the majority-spin 



(minority-spin) conduction electrons are in a metallic (in- 
sulating) state. We claim that this notion is not spe- 
cific to certain choices of system parameters but is fun- 
damental and ubiquitous for the ferromagnetic phase in 
the Kondo lattice model. Due to partial Kondo screen- 
ing, parts of the local moments are bound to the elec- 
trons, resulting in a dynamically-induced commensura- 
bility which is essential for producing the gap in the mi- 
nority spin electrons. We find that this commensurability 
condition leads to a nontrivial relation between electron 
magnetization, spin polarization and electron filling. 

The competition or cooperation between the magnetic 
phase mediated by the RKKY interaction and the Kondo 
screening can be modeled via a Kondo lattice model with 
antiferromagnetic coupling between the local moments 
and the conduction electrons. The Kondo lattice model 
reads [3, [TqI, El , 
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<i,j>a 



J y ^ SiSi 



where cj^ creates an electron on site i with spin-direction 
<j, p represents the vector of Pauli-matrices, and Si repre- 
sents the local spins which are coupled to the electrons via 
an antiferromagnetic spin-spin interaction with strength 
J > 0. 

To solve the Kondo lattice model we use the dynam- 
ical mean field theory (DMFT) DMFT maps 
the lattice model onto a quantum impurity model with 
a fermionic bath being determined self-consistently. Al- 
though being an approximation to real systems, DMFT 
has provided many insights into the physical properties 
and can even captures subtle differences in the lattice 
geometry. For solving the impurity model, we use the 
numerical renormalization group (NRG) 15, 16|, which 
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Figure 1: (Color online) Spin-resolved spectral functions for 
the ferromagnetic state in the Kondo lattice model J IW — 
0.25. The inset shows a magnification around the Fermi en- 
ergy for the spin-down component illustrating the gap in the 
spectral function. 



is able to reliably calculate spectral functions at very low 
temperatures [17|, |l8| . 

First, we briefly summarize the known DMFT results 

(A discussion 



n, m 2 



for the Kondo lattice model [K 
on the RKKY interaction within DMFT can be found in 
[l9|.) At half filling there is a pronounced antiferromag- 
netic Neel state for weak coupling, which vanishes with 
increasing coupling strength J via a continuous transition 
to a paramagnetic insulating state, the Kondo insula- 
tor. Doping slightly away from half filling this transition 
changes into a transition between an antiferromagnetic 
state (possibly spin-density-wave) and a paramagnetic 
metahic state. Especially the paramagnetic state around 
half filling is dominated by the Kondo effect, where the 
Kondo screening of localized spins results in a large Fermi 
surface accompanied by a narrow band and a gap close 
to the Fermi energy. Away from half filling the effects 
of Kondo screening become less important as there is an 
imbalance between local moments and available conduc- 
tion electrons. Such a tendency might be even stronger 
when the system enters a ferromagnetic state realized 
at low fillings, because an additional imbalance between 
spin- up and spin-down electrons arises. Contrary to this 
naive expectation, however, we demonstrate here that 
the Kondo screening plays an essential role even in the 
ferromagnetic phase. In particular, we reveal that the 
cooperation of ferromagnetism and Kondo screening can 
realize a novel kind of Kondo insulating state in the fer- 
romagnetic metallic phase. 

Figure [T] shows the local spin-resolved spectral- 
functions calculated in the ferromagnetic phase for a 
Bethe lattice with antiferromagnetic Kondo coupling 
JjW = 0.25 (bandwidth W = 4t). For this coupling 
strength the ferromagnetic phase extends from a nearly 
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Figure 2: (Color online) Gap width A/W in the minority- 
spin spectral function depending on the spin-coupling J and 
the occupation nj. The temperature of the system is T/W — 
3 • 10~^. The lines in the left panel are fits as ^ exp(— a/J). 



empty system to approximately = + = 0.5. 
One finds a striking difference in the spectral functions, 
which has not been recognized previously, for the major- 
ity spin (n|) and the minority spin (n|). While in the 
majority-spin spectral function a peak at the Fermi en- 
ergy u = and a dip for a; > can be found, there is 
a gap at the Fermi energy in the minority- spin spectral 
function. It is important to note that such a gap is not 
present in the ferromagnetic phase for a Kondo lattice 
model with ferromagnetically coupled spins. We propose 
that this gap in the spectral function is due to a partial 
Kondo screening of the localized spins, which results in 
an intriguing state: although the ferromagnetic state is 
metallic, only the majority-spin electrons contribute to 
the low-temperature properties, in particular transport. 
The minority spins, even though not completely depleted, 
form an insulator, which we name spin- selective Kondo 
insulator. 

Increasing the occupation number, the dip in the 
majority-spin spectral function moves closer to the Fermi 
energy and becomes more pronounced. Eventually, the 
ferromagnetic state is replaced by a paramagnetic state, 
for which the spectral functions for both spin-components 
suffer from the typical suppression of the DOS for > 
due to Kondo screening. Increasing the occupation to- 
wards half filling this dip becomes deeper and finally 
moves to the Fermi energy, forming the Kondo insula- 
tor. 

Clear evidence showing that the above insulating gap 
is indeed caused by the Kondo screening can be found 
in the dependence of the gap width A on the occupation 
and coupling strength, shown in Fig. O The left panel 
in Fig. [2] displays the dependence of the gap width on 
the coupling strength. For this purpose the minority- 
spin occupancy was kept constant (also resulting in a 
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Figure 3: (Color online) Upper panel: Magnetization and 
"commensurability", (nj + n;J), for two different coupling 
strengths and different occupation numbers calculated for a 
Bethe lattice at T/W — 3 • 10~^. The electron magnetization 
is shown as (m) = — nj, while the spin expectation value 
{Sz) is shown mirrored as —{Sz). The commensurability- 
condition is explained in the text. Lower panel: sketch of the 
local configuration (see text). 



nearly constant majority-spin occupancy). The depen- 
dency on J perfectly obeys a Kondo temperature-like 
form A ~ exp(— a/J) with a fitting constant a, sug- 
gesting that the Kondo physics is essential for the gap- 
formation. As a function of increasing filling the gap 
width decreases monotonically, as shown in the right 
panel of Fig. O As soon as the ferromagnetic phase van- 
ishes, the gap at the Fermi energy closes, too, and the 
minority and majority spectral functions look similar to 
the right panel in Fig. [H 

Let us now elucidate the basic physics behind this fer- 
romagnetic state. In Fig. [3] the magnetization of the 



conduction electrons (m) = n| — n| and the polariza- 
tion of the localized spins —{Sz) is shown. Note that 
the local spin-polarization always has the sign opposite 
to the conduction electron magnetization due to the an- 
tiferromagnetic coupling, thus —{Sz) has the same sign. 
Increasing the number of conduction electrons, the mag- 
netization of the electrons first increases due to increas- 



ing filling, and eventually decreases again due to the 
suppression of the ferromagnetic state. On the other 
hand, the spins are almost fully polarized for a nearly 
empty lattice, with monotonically decreasing polariza- 
tion for increasing conduction electron number. In the 
spirit of a pseudo-fermion representation, let us assume 
that the localized spins are actually formed by a local 
half-filled and strongly-interacting energy level so that 
{Sz) = {n^ — n;[)/2 and -\- n^^ = 1 (defining as the 
spin-dependent occupation of this level). Remarkably, 
we find that the following nontrivial commensurability 
condition holds within the ferromagnetic state: 

n| + nf = l, (1) 

as can be seen in Fig. [3l Note that Eq. ^ is equiva- 



lent to n| 



It should be noticed that this condition 



is not a priori given but is generated dynamically due 
to many-body effects. To clarify the origin of the above 
commensurability we propose that a partial local Kondo- 
singlet is formed in which (n|) majority- and minority- 
electrons participate, thus combining all spin-down con- 
duction electrons together with a part of the /-electrons 
and the spin-up conduction electrons to a Kondo spin- 
singlet. Here, we have assumed that spin-down is the 
minority-spin direction. The remaining majority-spin 
conduction electrons and spin-down /-electrons form a 
ferromagnetic state, (see a sketch in the lower panel of 
Fig. [3]). That the number of spin-down electrons in- 
cluding /- and conduction-electrons sums up to unity 
gives a commensurable situation, which results in a gap 
at the Fermi energy. On the other hand, for the majority- 
spins there is not such a commensurability condition, but 
-\- = = -\- n| holds. Therefore, this par- 
tial Kondo screening results in an insulating state for 
the minority-spin, while the majority-spin electrons re- 
main metallic. For this reason we have called this state 
a "spin-selective Kondo insulator" . The commensurabil- 
ity condition ([T]) smoothly connects to the Kondo insu- 
lator at half filling, suggesting that this ferromagnetic 
state should exist up to half filling. However, our results 
clearly show that there is a transition from this ferromag- 
netic phase to a paramagnetic state at electron fillings 
for > nferro- This is only possible, if the expectation 
value {Sz) jumps, leading to a discontinuous phase tran- 
sition at nferro • Our finding of the discontinuous transi- 
tion completely agrees with the recent analytical results 
showing that non-analytic terms prevent the continuous 
transition from a ferromagnet to a paramagnet [21 1. 

A further important consequence deduced directly 
from the commensurability condition ([T]) is a nontrivial 
relation between electron magnetization, spin polariza- 
tion and occupation number: 



2{Sz) + (m) = (n^) - 1. 



(2) 



This formula connects these three quantities which are 
otherwise independent from each other. By arranging it- 
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Figure 4: (Color online) Momentum- resolved spectral func- 
tions for the square lattice and J jW — 0.3, n — + nj = 
0.25, m — — nl = 0.1. The right side always shows a 
magnification of the left side around the Fermi energy cj = 
represented by the green line. From top to bottom the fig- 
ures show the majority-spin and the minority-spin spectral 
function, respectively. 




Figure 5: (Color online) Momentum-resolved occupation 
number n(k) (same parameters as in Fig. U). Left (right) 
panel shows the majority- (minority-) spin component. The 
dotted line represents the Fermi surface for non-interacting 
electrons (for the majority-spin this lines coincides with the 
shown surface). Note that for improving the contrast, the 
occupation for the minority-spin electrons is displayed in the 
interval n{k) G [0,0.3]. 



self in this way the system can gain an additional energy 
originating from the partial Kondo screening. Note that 
it should be possible to verify such a relation experimen- 
tally. 

The formation of the gap in the spectral function 
does not depend on the lattice geometry. For exam- 
ple, it can also be found in DMFT calculations for 
a two-dimensional square lattice. Figure [4] shows the 
momentum-resolved spectral functions for J/W = 0.3 
and n| + n| = 0.25 for a square lattice. The right 
panels are magnifications around the Fermi energy. In 
the minority-spin spectral function (bottom right) the 
gap can be clearly seen. While the majority-spin elec- 
trons are renormalized for > with a finite life- time, 
the minority-spin electrons are renormalized for < 0. 
This behavior is also visualized in Fig. [5l in which the 



momentum-resolved occupation number is shown. While 
for the majority-spin electrons the occupation number 
distribution looks like in the non-interacting case, for the 
minority-spin electrons this function is actually smeared 
out as compared to a reference non-interacting system, 
i.e. we indeed observe a large Fermi volume here. 

Using DMRG, we have confirmed that the spin- 
selective Kondo insulator together with the commensura- 
bility condition can be found for the ferromagnetic phase 
of the one-dimensional (ID) Kondo lattice model, too. In 
fact, previous calculations for the ferromagnetic ground 
state observed a magnetization 6*^0^ = 1/2(L— A^) [22] (L 
system length, Nc electron number), which supports the 
commensurability condition, and two separated bands in 
the spectral functions [23[. From these precise analyses 
of the ID model, we conclude that our finding is ubiqui- 
tous for the Kondo lattice model and not an artifact of 
DMFT. 

In conclusion we have clarified the physics behind the 
ferromagnetic metallic phase realized in the Kondo lat- 
tice model. We have demonstrated that the cooperation 
of ferromagnetism and partial Kondo screening results 
in an intriguing phase, here named spin-selective Kondo 
insulator, where an insulating state is stabilized for the 
minority-spin electrons while the majority-spin electrons 
are still metallic. We believe that the mechanism pro- 
posed here, the dynamically generated commensurabil- 
ity, should be generic for the ferromagnetic phase in the 
Kondo lattice models. It alternatively provides the non- 
trivial relation between the electron magnetization, spin 
polarization and occupation number, for which the sys- 
tem can gain a maximum of additional energy. The pro- 
posed relation between the macroscopic quantities might 
be confirmed in experiments. Good candidates in this 
context are ferromagnetic heavy fermion compounds, es- 
pecially compounds having a large Kondo temperature. 
For such compounds a verification of the above stated 
relation might be possible. Furthermore, spin-resolved 
transport measurements should show metallic majority- 
spin but insulating minority-spin electrons as well as a 
large Fermi surface for the minority-spin component. 
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